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An Isotropic Vector Field Decompoi^ition Method for Use in Scientific Compntations 


Backgromid of the Invention 


This invention relates generally to the field of scientific comptitation and simxdation, and 
more particularly to scientific computations which utilize vector field quantities. 

Computational methods of modeling and predicting physical phenomena have traditionally 
used a Cartesian or other three-tuple systems (collectively known as curvilinear coordinate 
systems) as a basis for describing the vector fields involved in modeling the said physical 
phenomena. The aforementioned three-tuple basis chosen for describing the vector fields then 
dictates the mathematical methodology used to express standard vector calculus operations. The 
algorithms to model the physical phenomena implement in a discrete manner the inherently 
continuous vector calculus operations as expressed in the associated curvilinear coordinate 
system. Current computational methods typically involve a finite-diflference approximation of the 
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partial derivatives involved in the evaluation of standard vector calculus operations. The 
aforementioned finite-difference approximation can be of varying order of accuracy (in a Talyor 
series expansion sense). Note also that while most computational methods colocate all the 
components of each vector field at every node within a grid, some computational methods stagger 
the components (to produce an uncolocated grid - for example the Yee grid, which is used in 
some computational electromagnetic methods). Furthermore, not all grids use a consistent amount 
of spacing between grid points, producing stretched and/or unstructured grids. 

The use of any standard set of curvilinear coordinates (and the associated three-tuple that 
expresses the vector field for that coordinate system) gives rise to grid anisotropicies. Said 
anisotropicies introduce directionally-oriented computational aberrations not associated with the 
underlying physical phenomena being studied. As an exanple, in Cartesian grid based 
computational electromagnetic solvers, the vector fields which represent the electric and magnetic 
fields propagate at different velocities along different trajectories within the grid - for instance, 
fields propagating along a Cartesian axis differ in velocity from fields propagatiag along a 
direction not aligned with a Cartesian axis. 


The primary object of the invention is to provide more accurate computerized simulations 
of physical phenomena by alleviating the inherent grid anisotropicy of Cartesian grid (and other 
standard grids) based simulations. 


Sunmiarv of the Invention 



Other objects and advantages of the present invention will become apparent from the 
following descriptions, taken in connection with the accompanying drawings, wherein, by way of 
illustration and example, an embodiment of the present invention is disclosed. 

An isotropic vector field decomposition method for use in scientific computations 
comprising the steps of: a computational grid modeled on a specific arrangement of nodes, such 
that each node is equidistant from it's twelve nearest neighbors (this grid is hereafter referred to as 
an isotropic vector matrix (IVM)), a vector field decomposition technique utilizing six vector 
components at every node within the associated IVM, and techniques for implementing standard 
vector calculus operations within the associated IVM. The preferred embodiment of the invention 
are algorithms for computational electromagentics, computational fluid dynamics, and other 
general scientific applications which involve vector field quantities* 

The drawings constitute a part of this specification and include exemplary embodiments to 
the invention, which may be embodied in various forms. It is to be understood that in some 
instances various aspects of the invention may be shown exaggerated or enlarged to facilitate an 
understanding of the invention. 


Brief Description of the Drawings 


Figure 1 is the basic VE cell with a chosen Cartesian orientation. 


Figure 2 is a two-frequency IVM. 
J: Figure 3 shows the six basis vectors of the VE cell. 

Figure 4 shows the a-plane of the VE cell. 


y Figure 5 shows the b-plane of the VE cell. 

S3 Figure 6 shows the c-plane of the VE cell 


Figure 7 shows the d-plane of the VE cell. 


Figure 8 shows the calculation of S3. 


Detailed Description of the Preferred Embodiments 


Detailed descriptions of the preferred embodiment are provided herein. It is to be 
understood, however, that the present invention may be embodied in various forms. Therefore, 
specijSc details disclosed herein are not to be interpreted as limiting, but rather as a basis for the 
claims and as a representative basis for teaching one skilled in the art to employ the present 
invention in virtually any appropriately detailed system, structure or manner. 

S Turning first to Figure 1 , there is shovra a basic cell whose 12 exterior vertices are all 

^2 equidistant fironi the center node. For the remainder of the description, the basic cell in Figure 1 is 
F called the VE cell (or just the VE). Txmiing now to Figure 2, there is shown a sanple version of 
m an isotropic vector matrix, which is built upon the basic VE cell. To acconq>lish an in^X)rtant 

function of this invention, there is shown in Figure 3 the six basis vectors of the VE cell: ej - eg. 

'^J For definitiveness, the sbc basis vectors in Figure 3 are now expressed in Cartesian coordinates (in 
O accordance with the Cartesian orientation shown in Figure 1): 

^3 
H 
^5 
^6 


= ax(A/3/2)-ay(l/2), 
= -ax(V3/2)-ay(l/2), 
= -ax(l / 2^3) - ay(l / 2) + az(V2 / V3), 
= ax(l/V3) + az(V2/V3), 
= -ax(l / 2^3) + ay(l / 2) + az(V2 / Vs), 
= ay. 


^ erf '1'^ 


Note that the basis vectors ei - eg have unit length as do the standard Cartesian unit 
vectors ax - a^. Also know that (from inspection): 


H = 



es = 


- 62 - «1, 

% = 

-ei 

-62. 


p;.; The basic VE cell has twelve vertices around the center node. We can use a six-tuple to 

' 2 reference all of these aforementioned twelve points, with each element of the tuple representing a 

M displacement in the corresponding unit vector direction* The center of the VE cell will be 

Si designated to have an address of (0,0,0,0,0,0), and the vertice designated with the address 

f : (1 ,0,0,0,0,0) represents the point at the tip of the basis vector ej , A shorthand notation for the 

% vertice at (1,0,0,0,0,0) is (+1). That is, we will use (-6) to represent the vertice at (0,0,0,0,0,-1). 

^0 We now describe a method to construct the vector calculus curl operator on the VE cell. 

For the development of this finite-difference style evaluation of the curl within the VE, we use 

2 two vector S and T, which are related by the equation S=VxT, which is to be evahiated at S(0), 
which requires values of the vector T at T(+l), T(-l), T(+2), T(-2), T(+3), T(-3), T(-h4), T(-4), 
T(+5), T(-5), T(+6), and T(-6). That is, we use the value of T at all twelve exterior points of the 
VE ceU to obtain an approximation of S at the center of the VE cell. Moreover, the vectors S and 
T will be vmtten is terms of the VE basis vectors; that is: 

S = 8161+8262 + 8363 + 8464 + 8565 + 8666, 
T = Tiei+T262 + T3e3 + T464 + T5e5 + T666. 


Within the framework of the VE cell, we want to evaluate S=VxT, which is defined (as 
shown in any advanced vector calculus text) as 

S=limdA->0(l/dA)[n^T.dI], 

where dA is the area enclosed by the contour of integration, and n is the unit vector is the 
direction which makes the right-hand side of the equation takes on its maximum value. The 
contour integral in the definition is commonly called the circulation; which is a term that will be 
used in the remainder of this description. To accomplish an important aspect of the invention, we 
will evaluate contours around the four hexagonal planes of the VE cell. The four hexagonal planes 
of the VE cell are labeled the a-plane, the b-plane, the c-plane, and the d-plane. The 
aforementioned contours on the four hexagonal planes of the VE cell include only the exterior 
points of the VE cell, with each plane containing a unique set of six points. The exterior vertices 
on each of the hexagonal planes that make up the four contours are (with the vertices in this 
order): 

a-plane: (+1) (-2) (+6) (-1) (+2) (-6) 

b.plane: (+1) (+4) (+5) (-1) (-4) (-5) 

c-plane: (-2) (+4) (+5) (-l)(-4)(-3) 

d-plane: (-5) (-3) (+6) (+5) (+3) (-6) 
For each of the four contours to be evaluated, note that the lefthand side of the curl involves three 
components of the S vector in the VE basis. We now evaluate the circulation around the a-plane 
that is shown in Figure 4, and take the dot product of the result against the unit vector normal to 
the a-plane, n^ (~^z ^ case). Thus 

(S3e3+S4e4+S5e5)^naKl/dA)[^T41]a, 

Note that 


e3-na=e4-na=e5-na=(V2/V3), 

thus 

(S3+S4+S5) = (V3/V2)(l/dA)[^.dI]a, 


= a', 


which defines the variable a'. That is, a' =<V3/V2)(l/dA)^T-dI]a. The evaluation of as it 

involves standard vector calculus techniques available in the literature, will not be discussed 
herein. 


The evaluation of the contour around the b-plane that is shown in Figure 5 results in the 
follovring equation: 

(S2+S3-S6)=b', 

where •dl]^. The evaluation of the contour aroxmd the c-plane that is shown 

in Figure 6 results in the follovraig equation: 

(Sl-S5-S6)=c^ 

where c'=(V3/V2)(l/dA)[^ •dlj^. The evaluation of the contour around the d-plane that is shown 
in Figure 7 results in the following equation: 

(Si-S2+S4)-d', 

where d'=(V3/V2)(l/dA)^T •dl](j. Evidently, as the equations from the four contours show, 
a'-b'+c'-d'=0. 
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In accordance with an important part of the invention, the values for the components of 
the resulting S vector can now be written: 

Si = (c^+d')/4, 
S2= (b'-d')/4, 

53 = (aW)/4, 

54 = (a^+dy4, 
85= (aV)/4, 
S6= (-bV)/4. 

which describes how to discretely approximate the vector calculus curl function on a VE cell, and 
within an associated isotropic vector matrix. In Figure 8 is shown the evaluation of the S3 
component of S - which involves the a' and W contours. The divergence of a vector field and the 
gradient can also be derived within the framework of this isotropic vector field decomposition 
method. 

While the invention has teen described in connection with a preferred embodiment, it is 
not intended to limit the scope of the invention to the particular form set forth, but on the 
contrary, it is intended to cover such alternatives, modifications, and equivalents as may be 
included within the spirit and scope of the invention as defined by the appended claims. 
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